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ABSTRACT 

The conceptual difficulties associated with a cosmological constant have led to the investi- 
gation of alternative models in which the equation of state parameter, w = p/p, of the dark 
energy evolves with time. We show that combining the supernova type la observations with 
the constraints from WMAP observations restricts large variation of p(z) at low redshifts. The 
combination of these two observational constraints is stronger than either one. The results are 
completely consistent with the cosmological constant as the source of dark energy. 

Key words: Cosmology: cosmic microwave background, cosmological parameters. 



1 INTRODUCTION 



the data from the high redshift super- 
became available, several independent con- 
indicated the existence of a cosmological constant 



Well before 
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straints 

I Bagla. Padmanabhan & Narlika rl Il99fj: lOstriker & Steinhardtl 
ll995t lEfstathiou. Sutherland & Maddoxl l 199(1) . In the last decade, 
observational evidence for an accelerating universe has be- 
come conclusive with almost all other possibilities bein g ruled 
out b y observations of h igh redshift super novae iTonrv et"ai] 
120031: iBarris et alJ 120031 : iRiess et alj |2004 and the cosmic 
microwave background radiation (CMBR) I Mel chiorri et alj 
2000; ISpergel et alJ [2003). The accelerated expansion of the 
universe requires either a cosmological consta nt or some form 
of dark e nergy iPadmanabhanl 120031: iPeebles & Ratr3 120031: 
Sahni & Starobinskv 2000; Ellis 2003) to drive the acceleration, 
with w = p/p < —1/3. Although a cosmological constant is 
the simplest solution from a phenomenological point of view 
(requiring just one fine tuned parameter), there is no natural 
explanation of the small observed value. This has led theorists 
to develop models in which a field, typically a s calar field, 
rovi d es the source of dark energy , e.g., quintessence jStei nhardt 



2003t iMacorra & Piccinefi] 12000: lUrena-Lopez & Matosl 1 200*01 
Gonzalez-DiazlteOOCtlRitis & Marinj200lHSen & SeshadriExiJ 



Rubano & Scudellaro 2002; Bludman & Roos 2002), k-essence 
I Armendariz-Picon. Mukhanov & Steinhardt 2001; Chiba 2002; 
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2003 


: IJassall 12003: 


Aguirregabiria & Lazkoz 2004; Sen 


2003 


: iGornietalJ 120031 iGibbons! 


20031 iKim. Kim & Kin 


2003 


; IShiu & Wassermanl 120021: 


Choudhurv et alJ l2002t 



Frolov, Kofman & Starobinsky 2002; Gibbons 2002), phan- 



tom fields JCaldwelj 120021 lHao & Lil l2003t iGibbonsI I2OO: 
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Noiiri & Odinstov 2003; Carroll. Hoffman & Trodden 2003; 
Singh. Sami & Dadhichl 120031: iFramptonl 120031: lHao & Lil 1 2003 : 
G(*nzi"le^Dia*|20(*^ ^ ; 
Cline. Jeon & MoordT 2003). branes lUzawa&Soda 2001; Jassal 
20031 lBurgesdl2003tl*vmtonll2003t lGonzalez-Diazll2000l) . etc. In 
the absence of significant spatial variation in the dark energy, the 
key difference between such models and the one with the cosmo- 
logical constant is that, in general, w is a function of redshift z in 
the former. Perturbations in dark energy also lead to observable 
signatures, though thes e can easily be confused w ith other p hysical 



Weller & Lewis 2003). There have also been a few proposals for 


unified dark matter and dark energy 




Padmanabhan & Choudhurv 


2002: ICardone. 


Troisi & Capozziello 


2004; 


Gonzalez-Diaz 2003) 



but these models are yet to be developed in sufficient detail to 
allow direct comparison with observations in a fruitful manner. If 
the current observations had excluded w = — 1 then one could 
have immediately ruled out cosmological constant as a candidate; 
but since this is not the case, direct exploration of w (z) at different 
redshifts, in order to check possible dependence of w on z, is of 
importance. Observations constrain the entire suite of parameters 
that describe cosmological parameters and while it is possible 
to choose other cosmological parameters so that ACDM is not 
allowed, these models have not been ruled out so far. 

It has been known for some time that supernova ob- 
servations and constraints from structure formation can be 
combi ned to put stringent limits on models for dark energy 
IPerlmutter. Turner & Whitel Il999h . Several attempts have been 
made in the past to constrain the equation of state for dark energy, 
along with other cosmological parameters, using the observations 
of galaxy clustering, temperature anisotro pies in the CMBR and 
the high redshift supernovae JPerlmutter. Turner & Whitel "19991 
Tegrrtark et alj 120031 rWang & Tegmarkl |2004 iHuterer & Coorayl 
20041: iMultamaki. Manera & Gaztanag al 120041: iMaccio et al J 
20031 iLinder & Jenkins! 120031: IPogosvan. Bond & Contaldill2003l: 
Lee & N j2003l:lLee. Lee and N J2003n~ ^ 
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This work, which is in the same spirit, focuses exclusively on 
constraining the variation of equation of state for dark energy by 
using a combination of such observations. We wish to constrain the 
variation of dark energy while keeping most of the other cosmolog- 
ical parameters fixed around their favoured values. In particular, we 
study the effect of a varying equation of state on the angular power 
spectrum of the CMBR fluctuations. Dark energy is not expected 
to be dynamically significant at the time of decoupling and — in 
fact — models in which t his is not true are pl agued by slow growth 
of density perturbations jB enabed & Bemar deaull200ll : lAmendolal 
l2003t IPeDeo. Caldwell & Steinhardtll2003l) . Nevertheless, evolv- 
ing dark energy will affect the features of temperature anisotropies 
in the CMBR in at least two ways: (1) The angular scale of features 
in temperature anisotropy, like the peaks, will change since the an- 
gular diameter distance depends on the form of w(z). (2) The inte- 
grated Sachs-Wolfe (ISW) effect will also depend on the nature of 
dark energy and its evolution, this effect is more relevant at small /. 
Thus observations of temperature anisotropies in the CMBR can be 
used to constrain the evolution of dark energy. Combined with the 
supernova observations, this allows us to put tight constraints on 
the equation of state of dark energy and its evolution. Our approach 
here is to use the full WMAP angular power spectrum in order to 
ensure that both the effects mentioned above are captured in the 
analysis. For lower multipoles (I ^ 20) signals from the ISW ef- 
fect, perturbations in dark energy and reionisation need to be disen- 
tangled, however the relative importance of this part of the angular 
power spectrum is limited as we use the full angular power spec- 
trum from WMAR Our aim is to demonstrate that the combination 
of WMAP observations and high redshift supernova observations 
is a very powerful constraint on variations in dark energy, certainly 
more powerful that either of the observations used in isolation. As 
far as we know, most previous attempts to constrain the dark en- 
ergy sector using WMAP data have not used the full angular power 
spectrum. 



2 VARYING DARK ENERGY 

Supernova observations, as well as the angular power spectrum of 
temperature anisotropies in the microwave background put geomet- 
ric constraints on dark energy. Thus these observations constrain 
the Hubble parameter H{£), which for the models under consider- 
ation here can be written as: 

H 2 (z) = H 2 [Q nr (1 + zf + Q DE {z)/ e ° E ] (1) 

where we have ignored the contribution of radiation to energy den- 
sity. Observed fluxes from supernovae and location of peaks in 
the CMB constrain distances to supernovae and to the surface of 
last scattering (Eisenstein & Wh itell2004l) . respectively. Dark en- 
ergy can be constrained only at redshifts where it contributes sig- 
nificantly to the energy density, therefore observations are sensitive 
to changes in q de (z) at low redshifts and variations at high red- 
shifts (z 3> 1) will be difficult to detect. Variation in g DE (z) can 
be written in terms of changes in the equation of state parameter 
w. Once we allow for variation of w, one is dealing with a function 
which, technically, has an infinite number of parameters. Given a 
finite number of observations, one can always fine tune such a func- 
tion. From a practical point of view, it is necessary (and often suffi- 
cient) to represent w(z) using a small number of parameters and do 
the analysis. The results wil l necessarily have some amount of pa- 
rameterization dependence l Bassett, Corasaniti & Kunz 2004) but 
this can be controlled by choosing different forms of parameteriza- 



tion and determining the range of variations in the results. In this 
work, we shall use the two parameterizations: 

w{z) = w +w 1 - ; p = l,2 (2) 
(1 + z)p 

For both p = 1,2 we have tu(0) = u>o;w'(0) = u?i but the 
high redshift behaviour of these functions are different: w(oo) = 
wo + wi for p = 1 while w(oo) = too for p = 2. Hence p — 2 
can model a dark energy component which has the same equation 
of state at the present epoch and at high redshifts, with rapid vari- 
ation at low z. Observations are not very sensitive to variations in 
w(z) for z> 1, hence return to the present value is not of critical 
importance. However, it does allow us to probe rapid variations at 
small redshifts. Forp = 1, we can trust the results only if wq + w 1 
is well below zero at the time of decoupling so that dark energy 
is not relevant for the physics of recombination of the evolution of 
perturbations up to that epoch. As we shall see, the allowed range 
of parameters does not contain such models. This constraint on the 
asymptotic value means that this parameterisation cannot be used 

to probe rapid variations at small redshifts. 

We will take Q de + O nr . = 1 jMelchiorri et alJ l2QO0t 
SDersel et al. I l2003h where density parameter for the total non- 
relativistic component is denoted by Q nr and the density param- 
eter for dark energy is fide- We will make assumptions regarding 
a few other parameters as well: (i) We assume that Qb = 0.05 
and (ii) h — 0.7. (iii) We also assume that the density fluctua- 
tions have the Harrison-Zel'dovich spectrum, i.e., n — 1. We do 
this since (a) there have been numerous studies in the past, restrict- 
ing the standard cosmological parameters and (b) we wish to fo- 
cus on the effect of wi 7^ 0. One can quarrel with these assump- 
tions but these are sufficient to give us an idea of the parameter 
space available for variation of dark energy. A complete study of 
these issues will need to address two questions: what is the param- 
eter space available for variation of dark energy if all parameters 
of interest in cosmological models are allowed to vary, this essen- 
tially addresses the issue of potential degeneracies of variation of 
dark energy with other parameters, and, if dark energy is allowed 
to vary then does that increase or change the allowed range of other 
param eters? We are carrying out a de tailed study to address these 
issues (Jassal, Bagla & Padmanabhan 2004). We will comment on 
dependence of our results on some of these parameters whenever 
relevant. 



For the supernova analysis, which is by now standard 



( Padmanabhan & Choudhury 2003; Choudhury & Padmanabhan 


2003 


Alametal. 2003, 2004; Wans 


etalJ 12004: 


Jonsson et al 


2004 


Jimenez 2003; Amendola & Ouercellini 2003; 


Jimenez et al 


2003 


Wandl2004 15alv & Diorgovskil 


2004; Gona 


20041 Lb assett 


12004 


Bassett, Corasaniti & Kunz 2004 


Corasaniti et al. 2004J), we 



combine the data on 230 supernovae published by To nry et al. and 
also for the 23 supernovae listed by Barris et al. iTonrv et alj2003t 
iBarris et all2003HRiess et all2004l) . To account for the uncertainty 
in measurement of cosmological redshift, we neglect the data points 
at redshift z < 0.01. The points with Av > 0.5 are neglected be- 
cause of the uncertainty in host galaxy extinction, after these cuts 
we are left with data for 194 supernovae. We compare the distance 
modulus predicted by theoretical models to the observed variation 
of apparent magnitudes of high redshift supernovae. We use the \ 2 
minimization method to find the best fit model. The supernova data 
favours dark energy to have an equation of state —1.3 ^ w(0) ^ 
—0.7 at the present epoch, with the lower values favoured for larger 
Q„ r - Significant variation in energy density of dark energy is al- 
lowed by supernova observations (see Figure |2j, e.g. the energy 
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density can change by more than a factor of 3 by z = 1. The re- 
sults we obtain are consistent with similar analysis in other papers 
I Padmanabhan_(^Chojjdhur3 |20 Q3l [c houdhurv & Padmanabhan 
2003 t lAlam et alj 120031 |2004 IWang et alJ |2004 Ijonsson etafl 
2004; Jimenez 2003; Amendola & Ouercellini 2003; Jimenez et al. 
2003 [ IWandl2004 iDalv & Diorgovskill2004 lGondhoo4 iBassettl 
2004 iBassett. Corasaniti & Kund |2004 ICorasaniti et alj l2004h. 
Use of supernovae sorted on the basis of da ta quality l Ton rv et alj 
120031: Barris et al. 2003; Riess et al. 2004) gives similar, though 
somewhat tighter constraints. 

Let us now add constraints from the WMAP observations. 
The location of peaks in the angular power spectrum of tempera- 
ture anisotropies (which depends on the angular diameter distance, 
which — in turn — depends on w(z)) provides the main constraint. 
Other effects like ISW provide a further constraint but their impact 
is limited as these are restricted to / ^ 20 for most models whereas 
we use the full WMAP angular power spectrum (I ^ 900). Given 
that the physics of recombination and evolution of perturbations 
does not change if we only modify w(z) within some safe limits, 
any change in the observed angular power spectrum will be due to 
a change in w(z). The location of peaks in the angular power spec- 
trum of temperature anisotropies will provide the main constraint 
and hence fix the effective value of w(z). We can define a w e ff 
by requiring that the angular size of a critical scale like the Hubble 
radius at the time of decoupling is same in the model with a given 
w(z) and for the model where w(z) — w e ff = constant. Ob- 
servations mainly constrain the value of w e ff. Thus if the present 
value wo < w e ff then we should have mi > 0, and similarly if 
wo > w e f f then mi < in the parameterisations we have used. 
Of course, if wi = then w e ff = wo. 

The priors we have considered here is Ho = 
70km s" 1 Mpc -1 , n = 1, Q B = 0.05, r = 0, 
-1.5 < wo sC -0.7, -2 < wi sC 2 and 0.2 sC fi M < 0.4. We do 
not vary many of these parameters as our main aim is to study the 
allowed variation in dark energy. If we allow variation of all pa- 
rameters at once, it will become difficult to understand the impact 
of individual variations. (We have varied some of the parameters 
individually, e.g., we have checked that allowing for reionisation 
with a prior r ^ 0.4 does not make a significant difference to the 
results, even though the best fit value of optical depth r is different 
for different models.) We are carrying out a detailed study where 
variation of most of these param eters is allowed and we also bring 
in more constraints Jassal. Bagla & Padmanabhan (2004). 

We do not take into account perturbations in dark energy in 
this paper. Perturbations in dark energy can modify our results to 
some extent. The effect of perturbtations depends on the detailed 
model of dark energy and cannot be incorporated in the model in- 
dependent aprroach we have used here. Including adiabatic pertur- 
bations in dark energy does not change our result significantly as 
the effect of these perturbations changes the angular power spec- 
trum only for I ^ 20. The effect of perturbations in dark energy 
becomes more relevant if we consider the allowed values of all cos- 
mological parameters in a varying dark energy model. 

To bring in these constraints from WMAP data, we need to 
compare the angular power spectrum expected in a given model 
with the observed spectrum iHinshaw et all2003h . We compute the 
theoretical angu lar power spec trum for each model using CMB- 
FAST ihttp://www.cmbfast.ord) . The comparison is carried out us- 
ing the likelihood program made available by the WMAP team 
IVerde et alj|2003l) . We normalize the angular power spectrum by 
maximizing likelihood in a comparison with the WMAP data. 




w(0) 




w(0) 



Figure 1. The top panel shows the allowed region at 99% confidence level 
in the w'(0) — u>(0) plane for Q nr = 0.3 for the parameterization in 
Eq.j2| with p = 1. The purple region is excluded by supernova observa- 
tions. WMAP constraints rule out the region in green at the same confi- 
dence level. The thick dotted line divides models that violate the strong en- 
ergy condition from those that do not. Models which preserve strong energy 
condition are on the top right of the dotted lines. The lower panel shows the 
same plot for p = 2. 

3 RESULTS 

Analysis of the high redshift supernova data shows that a large re- 
gion in the parameter space f2 nr — wo — wi is allowed at the 99% 
confidence level. For wi — 0, the favoured values of wo decreases 
(is more negative) with increasing fl nr . We will show detailed plots 
only for Q nr = 0.3 and then discuss the variation of allowed range 
of parameters with £l nr . We have studied large regions in param- 
eter space explicitly and have not relied on any expansion around 
the best fit model. 
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A fairly rapid evolution in the equation of state is allowed, 
so much so that w(z) ^ —1/2 at z ^ 0.5 is consistent with the 
supernova observations. Fig. Q shows the allowed region at 99% 
confidence level in the Wi — wa plane for fl nT = 0.3, regions in 
the parameter space that are ruled out are blanked out in purple. The 
99% confidence limit corresponds to a A\ 2 = 11.3 for the three 
parameters we have here. These results compare well with earlier 
studies that have addressed the issue of dark energy. We are not 
showing the full region allowed by supernova constraints as much 
of it is ruled out by WMAP observations. 

Adding WMAP constraints reduces the allowed region in pa- 
rameter space in a significant manner. The green region is ruled out 
by WMAP observations at 99% confidence level in Fig.Q We have 
A\ 2 = 13.28 for both p = 1 and p — 2. The orientation of the 
confidence level contours is consistent with the assertion that CMB 
observations essentially constrain iv e ff. Dotted lines within this re- 
gion mark the dividing line between models for which w(z) —1 
at all times (to the right and above the dotted line) and models 
which have w(z) < —1 at some redshift. These lines separate mod- 
els that violate the condition (p + p) ^ from those which do not. 
Equations for these lines are: 



w(0) = -1 (u)'(0) > 0) 

iw(0)+u/(0) = -l («/(0)<0) 



iu(0) = -1 
iu(0)+«/(0)/2 



(«f'(0) > 0) 

O'(o) < o) 



Observations do not discriminate between these two types of mod- 
els but we have indicated this line for the benefit of those who have 
a strong theoretical preference. 

In Fig. Q the top panel is for p = 1 and the lower panel is 
for p = 2. It is clear that the key features are the same in two pan- 
els. For both the parameterizations the region allowed by WMAP 
observations is much smaller than that allowed by the supernova 
observations (and the region occupied by non-phantom models is 
still smaller.) The lower bound from supernovae becomes stronger 
than the CMB constraint at small wo- The quantitative differences 
between the figures can be attributed to difference in w(z) at at a 
fixed z in the two parameterizations. 

To get a better feel, we have shown the allowed values of 
p DE (z) for dark energy as a function of redshift z for fl„ r = 0.3 
with just supernova data and by combining the constraints from 
WMAP in Fig|2| The purple region is ruled out by constraints from 
observations of high redshift supernovae alone. In cases where the 
allowed region extended beyond the range of parameters that we 
have studied, we used the largest allowed value within the range 
we had studied to get the envelope shown for this plot. For exam- 
ple, the allowed values of w'(0) for supernova constraints go be- 
yond the range that we have studied (—2 ^ tu'(0) ^ 2). This leads 
us to underestimate the region allowed by supernova constraints in 
fighter this does not affect the WMAP constraints. In this sense, 
our main conclusion — that WMAP observations put stronger con- 
straint on w(z) — is better than it appears from this figure. The 
allowed variation in p DE (z) is slower than that for matter, but con- 
siderable evolution is allowed by the supernova observations. The 
green region is excluded by WMAP constraints. For clarity, we 
have only marked those regions that are not already excluded by 
the supernova observations. Models that never have w(z) < — 1 
live between the broken lines. Models that are above the upper bro- 




Figure 2. The top panel shows those values of p(z) for 2 ^ z that are 
disallowed by the supernova and WMAP constraints for Q nr = 0.3. The 
purple region is ruled out by supernova observations and WMAP constraints 
limit the allowed region by also ruling out the green region. The remaining 
region contains allowed models that are parameterized by Eq.|2| with p = 
1. The allowed region for models that do not violate the condition (p+p) ^ 
is smaller and lies between the dotted lines. The lower panel shows the 
same plot for p = 2. 



ken line aer typically those that have w(0) < —1 but a positive 
w'(Q). Of course, a model that lives in the allowed region but has a 
different functional form for w(z) from that in Eq. J2j may not be 
allowed. The results for two parameterizations are very similar for 
z < 2. 
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Figure 3. The figure shows confidence levels in Q nr — wo plane for su- 
pernova and WMAP constraints (for both parameterisations) with w'(z = 
0) = 0. 99% confidence levels are plotted for both observational con- 
straints. The bold solid contour is 99% confidence level for WMAP. It is 
clear that the region allowed by the combination of two observational con- 
straints is much smaller than the one allowed by supernova data. 



4 CONCLUSIONS 

The observations we have considered constrain distances and the 
main constraint from CMB observations is on a reduced quantity 
w e ff which is the integrated effective value for the angular diam- 
eter distance. If w(0) < w e // (™(0) > w e ff) then w'(0) > 
(w'(0) < 0) is expected, indeed this explains the orientation of 
the allowed region in figure Q There is a strong degeneracy be- 
tween W e ff and fl nr , dw e f f /dflnr < (see Figure[3}- The cos- 
mological constant itself is ruled out at 99% confidence level for 
n nr > 0.375 and only phantom models survive beyond this, thus 
if other observational constraints were to rule out fl nr ^ 0.375 
then we will be forced to work with phantom models. (Note that 
the Cosmological constant and other non-phantom models are not 
ruled out by supernova observations for these values of fl nr ) But 
with current observations, such a constraint does not exist and the 
cosmological constant model is allowed. Hence observations do not 
require varying dark energy or phantom models even though such 
models are consistent with observations. The detailed dependence 
on Q nr and other parameters will be explored in detail in a later 
work Ijassal. Bagla & P admanab hanl2004|) . 

In our analysis we have not taken into account any pertur- 
bations in dark energy. This effect becomes increasingly impor- 
tant as w approaches zero, leading to suppre s sion of pe rturba- 
tions in matter iBenabed & Bemardeaul 1200 ll : I Amendolj Eool 
IPeDeo. Caldwell & Steinhardtll2003l) " This effect is relevant only 
at I ^ 20 in most models and hence its relative importance is 
small as we are using the full WMAP angular power spectrum here. 
Perturbations in dark energy also have a non-adiabatic component, 
but this cannot be incorporated in our analysis as it requires a de- 
tailed model for dark energy whereas we are working with param- 
eterised variations. We are carrying out a more detailed analysis 



where we allow parameters like n, h, £2b, t, etc. to vary. Prelim- 
inary results suggest that after marginalizing over other parame- 
ters the region in the w(0) — w'(0) allowed by CMB observa- 
tions increases, mainly on the lower side of the region shown in 
figureQ Supernova observations then provide stronger constraints 
for smaller w'(0) while the CMB observations continue to pro- 
vide much stronger constraints on large w'(0). Adding other obser- 
vational constraints from structure formation reduces the allowed 
region significantly as extreme variations of parameters allowed 
by WMAP data predi ct too much or too little structure formation 
( Jassal . Bagla & P admanabhan 2004). These observations are not 
very relevant in the present work as we are keeping most parame- 
ters fixed. 

Our analysis of models and observations shows that the al- 
lowed variation of p DE (z) is strongly constrained by a combina- 
tion of WMAP data and observations of high redshift supernovae. 
Indeed, given the allowed window as seen in Fig.|2| the cosmolog- 
ical constant seems to be the most attractive dark energy candidate. 
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